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We report new shell-model calculations of the isospin-symmetry-breaking correction, 5c, to su- 
perallowed 0"'' 0^ nuclear (3 decay. The most important improvement is the inclusion of core 
orbitals, which are demonstrated to have a significant impact on the mismatch in the radial wave 
functions of the parent and daughter states. We determine which core orbitals are important to 
include from an examination of measured spectroscopic factors in single- nucleon pick-up reactions. 
In addition, where new sets of effective interactions have become available since our last calcula- 
tion, we now include them; this leads to small changes in 5ns as well. We also examine the new 
radiative-correction calculation by Marciano and Sirlin and, by a simple reorganization, show that 
it is possible to preserve the conventional separation into a nucleus-independent "inner" radiative 
term, AX, and a nucleus-dependent "outer" term, S'^ We tabulate the new values for 5c, 5ns and 
5'ji for twenty superallowed transitions, including the thirteen currently well-studied cases. With 
these new correction terms the corrected J-t values for the thirteen cases are statistically consistent 
with one another and the anomalousness of the result disappears. These new calculations lead 
to a lower average Tt value and a higher value for Vud- The sum of squares of the top-row elements 
of the CKM matrix now agrees exactly with unitarity. 

PACS numbers: 23.40.Bw, 23.40.Hc 



I. INTRODUCTION 

Superallowed 0"*" — > 0"*" nuclear (3 decay currently pro- 
vides the most precise value for Vud, the up-down cl- 
ement of the Cabibbo-Kobayashi-Maskawa (CKM) ma- 
trix P, [3, 3 . This element is the key ingredient of the 
most demanding available test of CKM-matrix unitarity, 
a fundamental requirement of the electroweak standard 
model. To extract Vud from the experimental data, small 
theoretical corrections - of order ^1% - must be applied 
to take account of unobserved radiative effects as well as 
the isospin symmetry-breaking that occurs between the 
analog parent and daughter states of each superallowed 
transition 0, . Even though these corrections are very 
small, experimental measurements have by now reached 
such high precision that the uncertainty on Vud (±0.03%) 
is currently dominated not by experiment but by the un- 
certainty on these theoretical corrections. 

In the determination of Vud, an important strength 
of the nuclear measurements is that there are many 
O"*" 0+ transitions available for study, and currently 
there are thirteen of them, ranging from to ^^Rb, 
that have been measured with high precision. With so 
many, it becomes possible to validate the analysis pro- 
cedure by checking that all transitions individually yield 
statistically consistent results for Vud- Since the isospin- 
symmetry-breaking corrections depend on nuclear struc- 
ture, they differ from transition to transition and are 
particularly sensitive to this consistency test. Thus the 
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appearance of an anomalous result from any transition 
could signal a problem with the structure-dependent cor- 
rection for that case, a problem which might have impli- 
cations for other cases as well. 

In the most recent survey of superallowed 0+ 0+ 
transitions, which appeared in 2005 the results for all 
precisely measured cases - there were twelve at that time 
- were statistically consistent with one another. Today, 
there are thirteen such cases and they still form a sta- 
tistically consistent ensemble overall. However, recent 
precise Penning-trap measurements @, 01 of the Qec 
value for the superallowed decay of have left the 
result for that transition more than two standard de- 
viations away from the average of all other well-known 
transitions. This possible anomaly led us initially to re- 
examine the isospin-symmetry-breaking corrections for 
the ^^V transition, but what we learned from that reex- 
amination prompted us to a more general reevaluation of 
the corrections for other transitions as well. 

Our previous shell-model calculations for "^^V consid- 
ered six valence nucleons occupying the p/-shell orbitals 
outside a '^"Ca closed shell. This model space generated 
reasonable energies and spins for the known states in 
''^Ti, the daughter of ^^V. However, an important part 
of the charge-dependent correction depends on the ra- 
dial mismatch between the decaying proton in the par- 
ent nucleus and the resulting neutron in the daughter 
nucleus; but both these nucleons are bound to 45Ti. so 
the structure of that nucleus turns out to be important 
too. What is most striking about ^^Ti is that it has 
a 3/2+ state at an excitation energy of only 330 keV, 
which is strongly populated in single-nucleon pick-up re- 
actions like {p,d) and (■^He,Q;). Such low- lying sd-shell 
states can contribute to the structural parentage of the 
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initial and final states of the superallowed transition and 
consequently must affect the radial mismatch between 
them. This indicated to us that a complete calculation 
of the isospin-symmetry-breaking correction for the de- 
cay of ■^^V should include contributions from shells deeper 
than the pf shell. 

Two questions then arose. How many deeper shells 
need to be included and, if this effect is important for "'^V 
decay, how many other transitions will be similarly af- 
fected? In section|nT]of this paper, we address these ques- 
tions and settle on criteria for including deeper shells. 
Using these criteria - and incorporating more recent ef- 
fective interactions that have become available since our 
last work - we then re-evaluate the isospin-symmetry- 
breaking corrections for all transitions of relevance to 
the study of superallowed 0+ O"*"/? decay. For the 
cases with A < 38 the changes in the corrections arc 
very small - typically 0.03% - but for the heavier nuclei 
the changes can be as large as 0.2%. Most significantly, 
with the new calculated corrections, the result for ''^V is 
no longer anomalous. 

In section IIVI we incorporate recent improvements 
made by Marciano and Sirlin Q to the calculation of the 
radiative corrections for superallowed decays and then in 
section |V] wc apply both types of corrections - isospin- 
symmetry-breaking and radiative - to the current exper- 
imental data for superallowed decays. The result for Vud 
is changed appreciably, although it is still within quoted 
uncertainties of its old value, and the CKM-unitarity sum 
is improved. 



II. SUPERALLOWED BETA DECAY 

Superallowed Fermi beta decay between 0"'" states de- 
pends uniquely on the vector part of the hadronic weak 
interaction. When it occurs between isospin T = 1 analog 
states, the conserved vector current (CVC) hypothesis in- 
dicates that the ft values should be the same irrespective 
of the nucleus, viz. 



ft 



K 



GI\Mf\ 



const, 



(1) 



where K/{ncf = 2T:^h\n2/ {rrieC^f = (8120.278 ± 
0.004) X 10~^° GeV^^'s; Gy is the vector couphng con- 
stant for semi-leptonic weak interactions; and Mp is the 
Fermi matrix element. The CVC hypothesis asserts that 
the vector coupling constant, Gv, is a true constant and 
not renormalised to another value in the nuclear medium. 

In practice, Eq. ^ has to be amended slightly. Firstly, 
there are radiative corrections because, for example, the 
emitted electron may emit a bremsstrahlung photon that 
goes undetected in the experiment. Secondly, isospin is 
not an exact symmetry in nuclei so the nuclear matrix ele- 
ment, Mp, is slightly reduced from its ideal value, leading 
us to write: 



where Mq is the exact-symmetry value, which for T = 1 
states is Mq = v^- Thus, we define a "corrected" Tt 
value as 



Tt = ft{l + 5R){l-5c) 



K 



2G2(1 + AJ 



const, (3) 



where 5c is the isospin-symmetry-breaking correction, 8r 
is the transition-dependent part of the radiative correc- 
tion, and is the transition- independent part. For- 
tunately these corrections are all of order 1% but, even 
so, to maintain an accuracy criterion of 0.1% they must 
be calculated with an accuracy of 10% of their central 
value. This is a demanding request, especially for the 
nuclear-structure-dependent corrections. 

To separate out those terms that are dependent on 
nuclear structure from those that are not, we split the 
transition-dependent radiative correction into two terms. 



5'^ 
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of which the first, (5^, is a function only of the electron's 
energy and the charge of the daughter nucleus Z; it there- 
fore depends on the particular nuclear decay, but is in- 
dependent of nuclear structure. The second term, Sns, 
like Sc, depends in its evaluation on the details of nuclear 
structure. To emphasize the different sensitivities of the 
correction terms, we rewrite the expression for J^t as 



n ~ ft{l + 5'j,){l + S 



NS 



Sc) 



K 



2G2(1 + Av)' 



(5) 



(2) 



where the first correction in brackets is independent 
of nuclear structure, while the second incorporates the 
structure-dependent terms. 

From Eq. ([5]) it can be seen that a measurement of any 
one superallowed transition establishes a single value for 
Gv; moreover, measurements of many transitions pro- 
vides an excellent test of the validity of the whole analy- 
sis. Since CVC requires a unique value of Gv, all the ex- 
tracted values should be identical within experimental 
uncertainties. 

The /i-value that characterizes any /3-transition de- 
pends on three measured quantities: the total transi- 
tion energy, Qec'i the half-hfe, ti/2, of the parent state; 
and the branching ratio, i?, for the particular transi- 
tion of interest. The Q-EC-value is required to deter- 
mine the statistical rate function, /, while the half-life 
and branching ratio combine to yield the partial half-life, 
t. In 2005 we published a new survey of world data on 
superallowed 0+ O"*" beta decays All previously 
published measurements were included, even those that 
were based on outdated calibrations if enough informa- 
tion was provided that they could be corrected to modern 
standards. In all, more than 125 independent measure- 
ments of comparable precision, spanning four decades, 
made the cut. In the two years since the survey was 
closed another ten relevant pub lications have appeared 
@, 0, i, E [m, [13, [ll ElU and we have now incor- 
porated these results into our data base. Based on these 
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FIG. 1: Results from the 2005 survey [l[ updated with more recent published resuhs [1, 0, S [U, [H, [H, H [S O- 
The uncorrected ft values for the thirteen best known superallowed decays (left) are compared with the same results after 
corrections have been applied to obtain values. Here we have used the corrections calculated by us in 2002 [j, which were 
used in the original survey. The shaded horizontal band gives one standard deviation around the average n value. 



data for the thirteen most precisely known transitions, 
we obtain the ft values shown on the left side of Fig- 
ure [U then, by incorporating the corrections calculated 
by us in 2002 |1| and used in our 2005 survey [H , we ob- 
tain the corrected J-t values plotted on the right side of 
the figure. 

Obviously the calculated corrections do a remarkable 
job eliminating the considerable scatter that is evident in 
the /t- value plot on the left but is absent in the corrected 
!Ft values shown on the right. Overall, the statistical 
agreement among the J^t values is quite satisfactory, the 
normalized being 0.8. Thus, considering that the cor- 
rection terms were evaluated completely independently 
of these data, the consistency among the !Ft values can 
be taken as strong evidence that the correction terms are, 
in general, soundly based. 

However, there is a small but noticeable deviation from 
the average at ^^V (and possibly ^^Sc), which has only 
been revealed by the recent Penning-trap measurements 
of the transition Qec values. Though its statistical 
significance appears rather marginal in the figure, it must 
be remarked that the uncertainties quoted on these J^t 
values have been very conservatively determined. The 
measured data for each input parameter - (5£;c-value, 
half-life and branching ratio - were separately evaluated 
[H and, if the measurements were inconsistent with one 
another, the weighted-average uncertainty for that pa- 
rameter was increased to account for that inconsistency. 
In effect, for such cases, the original uncertainties quoted 
with the published measurements were all increased by a 
common "scale factor" that was large enough to restore 
statistical consistency among the measurements. (These 
scale factors are tabulated for each parameter in Ref . [l[ ; 



they range from 1 to 3.6.) This method, which is also 
used by the Particle Data Group , leads to final aver- 
age values that have a high confidence level but it does so 
at the cost of producing uncertainties that are in many 
cases larger than would result from a strict statistical 
average. 

With this method of analysis in mind, the excursion of 
the "^^V JFt value cannot be entirely ignored as a possible 
signal that the nuclear-structure-dependent corrections 
in this mass region are deficient. It certainly proved to 
be sufficiently provocative that we were led to the rceval- 
uation of correction terms that is reported here. 



III. ISOSPIN-SYMMETRY BREAKING 
CORRECTION, Sc 

For weak vector interactions in hadron states, the CVC 
hypothesis protects the decay amplitudes from strong- 
interaction corrections. However, there is a caveat. The 
CVC hypothesis also requires the hadron state to be an 
exact eigenstate of SU(2) symmetry (isospin). In nuclei, 
SU{2) is always broken, albeit weakly, by Coulomb in- 
teractions between protons. There may be other charge- 
dependent effects as well. These influences shift the value 
of the hadron matrix element from its exact symme- 
try limit to a new value and this shift has to be evalu- 
ated before weak-interaction physics can be probed with 
hadrons. In the case of superallowed (3 decay, the hadron 
matrix element, Mp, is given by Eq. ([2]) and it is Sc that 
we seek to evaluate. 

In the shell model for the cases of interest here, the A- 
particle wave functions representing the initial and final 
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states for supcrallowed f3 decay, |z) and |/), are states 
of angular momentum zero and isospin one. In a second 
quantisation formulation, the Fermi matrix element is 
written 

Mf = (/|r+|^) = Y,{f\aiap\z){a\T+\(3), (6) 

a, 13 

where the operator for Fermi f3 decay is the isospin ladder 
operator, aj^ creates a neutron in quantum state a and 
0/3 annihilates a proton in quantum state p. The single- 
particle matrix element, (a|r+|/3), is just a radial integral 



6^,p Rl{r)Rl{ry dr = 5^,p r^. (7) 



If the proton and neutron radial functions RJ^{r) and 
R^{r) arc identical, then the radial integral reduces to 
the normalization integral and has the value Tq, = 1. 

Now we introduce into Eq. ([6]) a complete set of states 
for the {A — l)-particle system, |7r), by writing 



(8) 



This is the essence of our model: we have allowed the 
radial integral to depend on the parentage expansion. 
Thus, we have added an additional label to and now 
write rj. 

If isospin is an exact symmetry, then the matrix ele- 
ments of the creation and annihilation operators are re- 
lated by hermiticity, (ttIoqIi) = (/|ajj,|7r)*. With that 
requirement, and with the radial integrals set to unity, 
the symmetry-limit matrix element is 



(9) 



Thus we see that the breakdown of isospin symmetry can 
enter the evaluation of Mp in one of two ways: either 
the matrix elements of and a^, are not related by her- 
miticity, or the radial integrals are not unity. Since each 
effect is small, we can, to first order, write the isospin- 
symmetry breaking correction as the sum of two terms 



(10) 



where in evaluating dci all radial integrals are set to unity 
but the matrix elements are not assumed to be related 
by hermiticity, while in evaluating Sc2 it is assumed that 
(TT\aa\i) = (/lall"")* but the radial integrals are allowed 
to diflFer from unity. Past calculations [1, [H] have indi- 
cated the radial overlap correction, dc2, is the larger of 
the two corrections so we will study this first. 



A. Radial Overlap Correction, 5c2 

1. Strategy for calculation 
For the Sc2 calculation, the Fermi matrix element is 



= Moh--L^|(/|al|7r)pf}M (11) 

where Mq is the exact-symmetry value, Eq. ([9|), and 
has been introduced as a radial-mismatch factor 



(12) 



Recalling that Sc2 is defined as |AfFp = |-^^oP(l — 5 02) 
we obtain 



(13) 



to first order in small quantities. A large contribution 
to (5(72 therefore requires a large spectroscopic amplitude 
and a significant departure of the radial integral from 
unity. 

There is an opportunity here to take guidance from ex- 
periment. The square of each spectroscopic amplitude, 
|(/|aj,|7r)p, is related to the spectroscopic factor mea- 
sured in neutron pick-up direct reactions. The exact 
relation, after inserting the isospin angular momentum 
couplings, is 



3C2 



E 



T^T; + 1) + 1 -r^(r^ + i) 

Tf{Tf + l) 



sly, (14) 



where .S'^'^j-^ is the spectroscopic factor for pick up of 
a neutron in quantum state a from an A-particle state 
of isospin T/ to an [A ~ l)-particle state of isospin Tt^. 
On setting Tf = 1 and separately identfying sums to 
the isospin- lesser states with T^. = |, denoted tt^, and 
the isospin-greater states with Tjr = |, denoted 7r>, we 
obtain a very revealing formula 



^C2 



E^^ 



E 



(15) 



This equation provides the key to the strategy we will 
use in calculating Sc2- It demonstrates that there is 
a cancellation between the contributions of the isospin- 
lesser states and the isospin-greater states. Moreover, if 
the orbital a were completely full in the initial A-particle 
wavefunction, then the Macfarlane and French sum rules 
[l9t for spectroscopic factors would require = 
1 '^TT> Sa ^^'^ the cancellation in Eq. ((T5|) would be very 
strong. In fact, the cancellation would be complete if 
fi^ ~ fl^. As wc will discuss further in the next section, 
this cancellation is not in general complete because the 
radial-mismatch factors for isospin-lesser states are larger 
than those for isospin-greater states. Even so, cancella- 
tion is always significant, and it becomes most complete 
when closed-shell orbitals are involved. Furthermore, the 
more deeply bound the closed-shell orbital, the greater 
the energy spread in the spectroscopic strength and the 



5 



TABLE I: Illustration of the strategy used in calculating Sc2 for V. The measured spectroscopic factors from the 
*®Ti('^He, a)*^Ti reaction [l3| are shown for the states where they are largest. Two calculations are then given for each 
state's contribution to 5c2' the first assumes that the total Macfarlane- French (M-F) sum rule is exhausted in each state, while 
the second gives the result of a complete shell-model calculation. Both methods give remarkably similar results. 
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more complete the cancellation. Thus, although the dom- 
inant contributions to Sc2 come from unfilled orbitals, we 
conclude that closed-shell orbitals must play a role, albeit 
one that decreases in importance as the orbitals become 
more deeply bound. 

Based on these observations, our strategy is to use ex- 
periment to guide us in determining which closed-shell or- 
bitals are important enough to include. Ideally, of course, 
one would take the spectroscopic factors determined from 
experiment and insert them into Eq. (jl5[) but, especially 
where delicate cancellations are involved, the reliability 
of (forty-year-old) experimental spectroscopic factors is 
certainly not up to the task. Our strategy then is to 
use the shell model to calculate the spectroscopic ampli- 
tudes in Eq. (|13p but to limit the sum over orbitals a just 
to those for which large spectroscopic factors have been 
observed in neutron pick-up reactions. 

We illustrate the strategy for the case of V. The spec- 
troscopic factors for neutron pick up from "^^Ti have been 
measured in the ('^He, a) reaction by Borlin [l^. He iden- 
tified sixteen states in ^^Ti, and in Table U we record 
the six states with the largest spectroscopic factors, i.e. 
S > 0.5. We note that the errors on the experimental 
spectroscopic factors are quite large, and in two cases 
the quoted Sa values (column 4) exceed the Macfarlane- 
French sum rule [loj for pure configurations (column 6). 
Thus we do not use the experimental spectroscopic fac- 
tor explicitly, but take them as a guide for which orbitals 
should be included in the shell-model calculation. In the 
case of '*^V decay, they tell us that orbitals /7/2, c?3/2 
and Sl/2 should be included. In column five of Table U 
we give a typical value for the radial mismatch factor, 
rij, for the given orbital a and isospin T^^. Column seven 
gives the contribution to Sc2 fi'om this a and isospin T^r 
if the Macfarlane- French sum rule is used for the spectro- 
scopic factor, while in columns eight and nine are shown 
the results of a detailed shell-model calculation. The re- 
sults from the Macfarlane- French sum rules and the shell- 



model calculation are remarkably similar. The summed 
Sc2 for the shell-model calculation (the sum of all en- 
tries in column 9) is 0.58%, nearly a factor of two larger 
than our previous calculated value, which was published 
in 2002 [|. 

The difference between our calculations arises as fol- 
lows: In 2002 our shell-model calculations for ^^V were 
based on the model space (fp)^, with six valence nucle- 
ons occupying the p/-shell orbitals. In fact, only the fj/2 
orbital contributed importantly to the Sc2 calculation so 
the result was 6c2 = 0.45 — 0.12 = 0.33% (see the two 
rows for the /7 /2 orbital in Table |T| . Absent from this 
2002 calculation was any contribution from the core or- 
bitals, Sl/2. In our present calculations, these 
orbitals are included, with the (1^/2 orbital contributing 
0.14 % to Sc2 and the si/2 contributing 0.11 %. 

But why stop there? Why not include the c?5/2 and 
possibly the p-shell orbitals in the computation? Our 
answer is that the neutron pick-up measurement saw lit- 
tle or no evidence for such core states, which implies 
that their spectroscopic strength is distributed widely 
over many states. In this case, the cancellation between 
isospin-lesser and isospin-greater states becomes more 
complete and their contribution to Sc2 is reduced to a 
level that wc believe can be neglected. 

With this approach, we are now in a position to re- 
vise our earlier results [J to include the effects of previ- 
ously ignored core orbitals. Again using measured spec- 
troscopic factors from neutron pick-up reactions, we de- 
termined that changes were required for the ^ = 22 and 
26 cases, in which p-shell holes must contribute in ad- 
dition to the original sc?-shell configurations; similarly, 
sc?-shell holes were required in addition to the p/-shell 
particles for A = 46, 50 and 54. For A = 62, 66, 70 
and 74 in the upper p/-shell there are no experimental 
neutron pick-up reaction measurements to guide us. Our 
previously published calculations for these nuclei were 
based on (P3/2, /5/2J Pr/2)" model spaces using ^^Ni as 
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a closed-shell core. It seemed prudent now for these cases 
at least to include the /7/2 orbital in the calculation of 
i5c2; and we have made this change. In the cases with 
A = 18 and 42, we had previously included some con- 
tribution from deeper shells; we did not need to make 
any changes in the former but did add the Si/2 and (^5/2 
shells to the latter. No additional orbitals were required 
for the cases with A = 10, 14, 30, 34 and 38. 



2. Radial-mismatch Factor, fi^ 

In considering the radial integrals, we benefit from a 
very strong constraint: the asymptotic forms of all radial 
functions must match the measured separation energies, 
Sp and Sn, where Sp is the proton separation energy 
in the decaying nucleus and S„ the neutron separation 
energy in the daughter nucleus. The basic ingredients 
of these separation energies arc well known and can be 
found in any atomic mass tables. It is the size of the dif- 
ference between Sp and Sn and the presence or absence of 
nodes in the radial wave functions that arc the principal 
factors in determining the magnitude of f2^. 

Our calculations of this mismatch factor follows the 
same path as that described in our earlier works [1, [l^l ■ 
We use a Saxon- Woods potential defined for a nucleus of 
mass A and charge Z + 1 as: 

V{r) = -Fo/(r) - V,g{r)\.a- + Vc{r) - Vgg{r) - V,Mr), 



where 



(16) 
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Ze^ ( 

3 , f„, ,.<«.. (17, 



with R = ro{A - 1)1/3 and R, = r,{A - l)i/3. The 
first three terms in Eq. are the central, spin-orbit 
and Coulomb terms respectively. The fourth and fifth 
terms arc additional surface terms whose role wc discuss 
shortly. 

Most of the parameters were fixed at standard values, 
Vs = 7 MeV, Ts = 1.1 fm and a = Us = 0.65 fm. The 
radius of the Coulomb potential was determined from the 
charge mean square radius, {r'^)lC^ , of the decaying nu- 
cleus as determined from elastic electron scattering; see 
Eqs. (21) and (22) in Ref. 0]. The well radius, Tq, was 
similarly fixed, by requiring that the charge density con- 
structed from the square of the proton wave functions 



bound in the well should also match the charge mean 
square radius. Initially, with Vg and Vh set to zero, the 
well depth, Vq, was adjusted so that the binding energy 
of the least-bound orbital matched the experimental sep- 
aration energy. 

Prom the shell model calculation, wc obtained the A- 
particle wave functions, \i) and \f), expanded into prod- 
ucts of {A — l)-particle wave functions |7r) and single- 
particle functions |a). In Eq. ([8]) and the discussion that 
followed it, we noted that the radial integral should de- 
pend on the separation energies relative to the (^ — 1) 
state, \tt). We ultimately allowed this to happen but ini- 
tially we calculated the value of Sc2 under the assumption 
that the proton and neutron radial functions, R^{r) and 
R"{r), have asymptotic forms for all a that are fixed at 
the separation energies, Sp and Sn, to the ground state 
of the {A — 1) nucleus. In this case, the sums over tt 
can be done analytically and the computed value of 5c2 
becomes independent of the shell-model effective inter- 
action. This result, which we label S^2, can be simply 
expressed with the help of Eqs. ^ and 



S'r. 



C2 



(18) 



Here ttg is the shell-model orbital of the transferred neu- 
tron in the pick-up reaction from the A-particle state |/) 
to the ground state of the {A — l)-particle nucleus. 

We next removed our simplifying assumption and eval- 
uated the radial integrals with eigenfunctions of the 
Saxon- Woods potential whose well depth was adjusted 
so that each eigenfunction matched the separation en- 
ergy of the {A — 1) state to which it corresponds, [tt). 
For an (A — 1) state at excitation energy the corre- 
sponding separation energies are Sp + and Sn + E^. 
We label these results 6^2 and note that the values now 
depend on the spectroscopic amplitudes, and hence on 
the shell-model effective interaction, but not strongly. 

So far, we have ignored the two surface terms in 
Eq. pT|) by setting = and Vh = 0. It can be argued, 
however, that the central part of the potential, which 
in principle should be determined from some Hartree- 
Fock procedure, should not be continually adjusted. In- 
stead, any adjustments made to match separation ener- 
gies should be to the surface part of the potential rather 
than to the depth of the well. Thus, we also calculated 
Sc2 by fixing Vq separately for protons and neutrons to 
match the ground-state parent separation energies, Sp 
and Sn, and then adjusting the strength of the surface 
term, Vg (keeping Vh = 0) so that the asymptotic forms 
matched the separation energies Sp + E^ and Sn + E^. 
These results are labelled 5{J2 ■ 

Finally, our fourth method of calculation was the same 
as the third, except that it was the second surface term, 
Vh, that was adjusted to match separation energies, keep- 
ing Vg = 0. This second term, h{r), is even more strongly 
peaked in the surface than g{r). These results are la- 
belled (5^^. 

On average, the method III values of Sc2 are about 2% 
lower than the method II values; and method IV values 
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are about 7% lower than the method II values for orbitals 
without any radial nodes. For orbitals with one or more 
nodes, there is more of the radial wave function in the 
surface region and methods III and IV produce greater 
reductions. 



TABLE II: Calculations of 5c2 with Saxon- Woods radial func- 
tions, without parentage expansions (5^2 ) ^-nd with parentage 
expansions ((5^21 ^C2 j s-nd Sq^). Note that only one sample 
result is shown in each case for 1 ^C2i ^c-I S'ld 5q^, while 
the adopted 5c2 value in column 7 reflects the results from 
all multiple-parentage calculations for that case; see text. 



3. The shell-model calculations 

We now present our results for 5c2 based on the ex- 
tensions of the shell-model spaces mentioned at the end 
of Sect, nil A II In addition to adding the core orbitals 
mentioned there, however, in some cases wc have also 
been able to make use of more recent effective interac- 
tions that have become available since our last work. 
Specifically, we have used the following interactions in 
the various mass regions of interest: In the p-shell, we use 
the Cohen-Kurath interactions [2l| and the more recent 
PWBT interaction of Warburton and Brown (2^ . In the 
s, d-shell, besides the universal interaction of Wildenthal 
[231, we employ two new versions, USD-A and USD-B, of 
Brown and Richter [23| . In the p/-shell we use the KB3 
interaction of Kuo-Brown [2^ as modified by Poves and 
Zuker [2^ , the FPMI3 interaction of Richter and Brown 
[131, and the more recent GXPFl interaction of Honma 
et al. [H, [2^ . For cross-shell interactions between the 
major shells, we have used the interaction of Millener 
and Kurath [3^. It should be noted that in many cases 
we found it necessary to introduce some truncations in 
the original model space in order to keep the calculations 
tractible. 

We made calculations for all twenty superallowed tran- 
sitions considered in our earlier work and for each 
we calculated 5c2 in the four methods, I-IV, described 
in Sect. lIII A 21 and with the several interactions listed 
in the previous paragraph. In Table [U we record only 
one sample result for (5^2: <^C2' ^ci ^"^^ ^C2 ^'-'^ each 
nucleus listed. However, our "adopted 602" values re- 
sult from our assessment of all multiple-parentage cal- 
culations made for each decay, not just those shown in 
the previous three columns. The uncertainty assigned to 
each adopted value reflects the uncertainty in the radius 
of the Saxon- Woods potential (resulting from an uncer- 
tainty in the nuclear rms radius to which it is adjusted), 
the spread of results obtained with different shell-model 
interactions, and the spread of results obtained with the 
different procedures labelled II, III and IV in the Table. 



B. Isospin-Mixing Correction, Sci 

The second (and smaller) contribution to 5c is the 
isospin-mixing correction, dci- For its evaluation, the 
radial integrals are all set to unity, but the spectroscopic 
amplitudes in Eq. ([8]) are not required to satisfy hermitic- 
ity. Calculations of this correction turn out to be very 
sensitive to the details of the shell-model computation. 
This would be a very unfortunate property if wc were 





2002 


This work 




Parent 
nucleus 


5C2{%) 

Ref. \£ 


S'c2i%) -502 (%) 5c"(%) 5C2(%) 


5c2i%) 

adopted 



10 c 

14 Q 

i«Ne 
22Mg 
2«Si 

30 g 

Ar 

38 Ca 

42 



-1 : 

0.170(15) 
0.270(15) 
0.390(10) 
0.255(10) 
0.330(10) 
0.740(20) 
0.610(40) 
0.710(50) 
0.555(40) 



0.132 
0.217 
0.251 
0.207 
0.223 
0.812 
0.351 
0.402 
0.359 



0.163 
0.274 
0.386 
0.366 
0.421 
0.714 
0.680 
0.840 
0.881 



0.165 
0.271 
0.387 
0.382 
0.407 
0.710 
0.639 
0.784 
0.849 



0.163 0.165(15) 

0.271 0.275(15) 

0.382 0.410(25) 

0.375 0.370(20) 

0.392 0.405(25) 

0.713 0.700(20) 

0.579 0.635(55) 

0.702 0.745(70) 

0.780 0.835(75) 



= 



26 ^1"^^ 


0.230(10) 


0.156 


0.292 


0.280 


0.271 


0.280(15) 


3*01 


0.530(30) 


0.312 


0.583 


0.561 


0.498 


0.550(45) 


38 j^m 


0.520(40) 


0.299 


0.623 


0.575 


0.522 


0.550(55) 


^^Sc 


0.430(30) 


0.278 


0.681 


0.648 


0.606 


0.645(55) 


46y 


0.330(25) 


0.273 


0.587 


0.543 


0.506 


0.545(55) 


^^^Mn 


0.450(30) 


0.315 


0.638 


0.598 


0.594 


0.610(50) 


s-^Co 


0.570(40) 


0.376 


0.760 


0.688 


0.706 


0.720(60) 


«2Ga 


1.05(15) 


1.31 


1.22 


1.19 


1.14 


1.20(20) 


"^As 


1.15(15) 


1.32 


1.41 


1.34 


1.24 


1.35(40) 


™Br 


1.00(20) 


1.43 


1.41 


1.31 


1.10 


1.25(25) 




1.30(40) 


1.68 


1.60 


1.47 


1.12 


1.50(30) 



not able to adopt certain strategies that act to reduce 
the model dependence considerably. 

There are three ways in which we incorporated charge 
dependence in our shell-model calculation. First, the 
single-particle energies of the proton orbits were shifted 
relative to those of the neutrons. The amount of shift was 
determined from the spectrum of single-particle states in 
the closed-shell-plus-proton versus the closed-shell-plus- 
neutron nucleus, where the closed shell was taken to be 
the nucleus used as a closed-shell core in the shell-model 
calculation. Wc took these single-particle shifts from ex- 
periment and did not adjust them. 

Second, we added a two-body Coulomb interaction 
among the valence protons and adjusted its strength so 
that the measured 6-coefficient of the isobaric multiplet 
mass equation (IMME) was exactly reproduced. Third, 
we introduced a charge-dependent nuclear interaction by 
increasing all the T = 1 proton-neutron matrix elements 
by about 2% relative to the neutron-neutron matrix el- 
ements. The precise amount of this increment was de- 
termined by requiring agreement with the measured c- 
cocfficient of the IMME. This strategy of constraining 
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TABLE III: Shell-model calculations of the isospin-mixing correction, Sci- 





Measured 


2002 








This work 




Parent 


IMME coefficients [33.1 


fci(%) 


E40+) 




■) 5ci{%) 


5ci{%) 


5ci{%) 


nucleus 


b (keV) c (keV) 


Ref. \4} 


expt 


SM 


unsealed 


scaled 


adopted 



n = -1: 





-1.546 


0.362 


0.010(10) 


6.18 


9.24 


0.005 


0.011 


0.010(10) 


14 Q 


-2.493 


0.337 


0.050(20) 


6.59 


6.64 


0.049 


0.050 


0.055(20) 


"Ne 


-3.045(1) 


0.347(1) 


0.230(30) 


3.71 


4.07 


0.116 


0.140 


0.155(30) 


22 Mg 


-3.814(1) 


0.315(1) 


0.010(10) 


6.24 


6.21 


0.010 


0.010 


0.010(10) 


26Si 


-4.535(2) 


0.302(2) 


0.040(10) 


3.59 


3.86 


0.022 


0.026 


0.030(10) 


30 g 


-5.185(2) 


0.275(2) 


0.195(30) 


3.79 


3.80 


0.137 


0.138 


0.155(20) 


34 Ar 


-5.777(2) 


0.286(2) 


0.030(10) 


3.92 


3.97 


0.023 


0.023 


0.030(10) 


38Ca 


-6.328(3) 


0.284(3) 


0.020(10) 


3.38 


3.21 


0.026 


0.023 


0.020(10) 




-6.712(3) 


0.287(3) 


0.220(100) 


1.84 


3.16 


0.038 


0.114 


0.100(20) 



T, = 0: 





-4.535(2) 0.302(2) 


0.040(10) 


3.59 


3.86 


0.025 


0.028 


0.030(10) 


3*01 


-5.777(2) 0.286(2) 


0.105(20) 


3.92 


3.97 


0.091 


0.093 


0.100(10) 


38 j^m 


-6.328(3) 0.284(3) 


0.100(20) 


3.38 


3.21 


0.099 


0.089 


0.105(20) 


42Sc 


-6.712(3) 0.287(3) 


0.060(30) 


3.30" 


5.05 


0.007 


0.017 


0.020(10) 


46y 


-7.327(10) 0.276(11) 


0.095(20) 


3.57" 


4.86 


0.040 


0.075 


0.075(30) 




-7.892(30) 0.259(30) 


0.055(20) 


3.69 


3.62 


0.057 


0.054 


0.045(20) 


54Co 


-8.519(25) 0.276(25) 
-9.463(70) 0.265(25)'' 
-9.95(15) 0.262(25)'' 
-10.48(23) 0.260(25)'' 
-10.82(25) 0.258(25)'' 


0.040(15) 


2.56 


2.26 


0.058 


0.045 


0.050(30) 


62Ga 


0.330(40) 


2.33 


2.32 


0.221 


0.219 


0.275(55) 


66 As 


0.250(40) 


2.17^= 


1.89 


0.210 


0.159 


0.205(45) 


^°Br 


0.350(40) 


2.01 


2.05 


0.332 


0.346 


0.350(40) 
0.130(60)^* 


^4Rb 


0.130(60) 


0.508 


0.523 


0.122 


0.129 



" Second excited 0+ state; shell-model calculations indicate this state talics up most of the depletion from the analog state. 

^ Estimated: extrapolated from a fit to c coefficients in 0+ states in A = 4n -|- 2 nuclei, 10 < j4 < 58; the data were taken from Ref. [ssll . 

Estimated: value is the average of the excitation energy of the 0+ states in 62 Zn and ™Se. 

No new calculations were performed for ^'^Hh. 



the charge-dependence in the effective interaction by re- 
quiring it to reproduce the coefficients of the IMME was 
adopted from the work of Ormand and Brown [Sll [3^. 

Experimental data were used in one more way to con- 
strain our calculations. If isospin were an exact symme- 
try, then the parent 0+ (T = 1) state would decay exclu- 
sively to its analog state in the daughter nucleus. Beta 
transitions to all other 0+ states in the daughter would 
be strictly forbidden. But, with isospin symmetry bro- 
ken, weak transitions (with branching ratios measured in 
parts per million) can occur to these other 0+ states. In 
this case, we write the Fermi matrix element squared to 
the n*'' non-analog 0''' state as 



9A" 



(19) 



and the reduction in the analog transition Fermi matrix 
element squared as 



(20) 



neglecting, in this context, the contribution of Sc2- If ah 
the 0"'' states of a given model space had the same T = 1 



isospin designation, then the effect of isospin-symmetry 
breaking terms in the Hamiltonian would be to deplete 
the analog-transition strength by an amount that is ex- 
actly matched by the sum of the strengths to the non- 
analog states: i.e. 



3C1 



rn 

"CI- 



(21) 



In practice, with large shell-model calculations the 0+ 
states in the model space will include some states whose 
isospin designation is not T = 1 ; and Eq. (pij) is not then 
exactly correct. Nevertheless, it remains approximately 
true. 

Significantly, in many cases the bulk of the analog state 
depletion shows up in a single excited 0"*" state, usually 
(but not always) the first excited one. This allows us 
once again to use experiment to constrain and refine our 
calculation. In the limit of only two-state mixing, per- 
turbation theory would indicate that 
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where AE' is the energy separation of the analog and 
non-analog 0"*" states. Again, this is not an exact result, 
but it does highlight the importance of the shell-model 
Hamiltonian producing a good quality spectrum of 0^ 
states with, in particular, the first excited non-analog 
0+ state calculated to have an excitation energy close to 
its experimental value^. This is not always possible to 
achieve in the shell model, especially near closed shells 
where excited 0+ states tend to exhibit strong deforma- 
tions. We used two strategies to bring the calculation 
into line with experimental information. Our first was 
to adjust the centroids of the shell-model Hamiltonian 
matrix elements specifically to get the excited 0+ state 
at about the right energy. Our second was to scale our 
calculated Sci value by a factor {AE)^^^^^/ {AE)l^p^, the 
ratio of the square of the excitation energy of the first 
excited 0+ state in the model calculation to that known 
experimentally. 

We list in Table Hill the experimental values [1^ of the 
IMME coefficients, b and c, and the known excitation en- 
ergy Ex{0^) of the first (or second) excited 0+ state in 
the daughter nuclei. As explained, all our shell-model cal- 
culations were adjusted to reproduce exactly the values 
of b and c, and to match, as closely as possible the exci- 
tation energy of the excited 0+ state. We compensated 
for any remaining discrepancies between the calculated 
and experimental values of Ex{0^) by scaling the results 
for Sci- As in Table HIl we give (in columns 6-8) the re- 
sults from one sample calculation for each nucleus. Then 
in column nine we present adopted Sci values that result 
from our assessment of the results of all calculations made 
for each decay, not just the ones shown in columns 6-8; 
the uncertainties were chosen to encompass the spread 
in the results from those calculations and to include the 
uncertainty in the IMME b and c coefficients. For com- 
parisonjin column 4 we list the values we adopted for Sci 
in 2002 Our strategies have remained unchanged, but 
here we have additionally used some more recent shell- 
model effective interactions as listed in Sect. IIII A 31 In 
nearly all cases, the new values of 6ci agree with the old 
values within their stated uncertainties. 

For the heavier nuclei there are experimental data on 
Fermi transitions to the non-analog excited 0+ states. 
The measured branching ratios 15, 34, 35,, ,36| have been 
converted to S^^i values, via Eq. (jl9[) . and listed in Ta- 
ble llVI Again, for each nucleus, we list just one represen- 
tative calculation and our adopted value. The assigned 
error reflects both the spread among the different calcu- 
lations and the uncertainties in the IMME coefRcients. 
Our 2002 adopted values Q are also listed. For nuclei 
38 < A < 54, with the possible exception of ^°Mn, the 
agreement between theory and experiment is entirely sat- 



TABLE IV: Shell-model calculations of Sq^ for Fermi decay 
to the first excited 0^ state; see Eq. (|19p . The results are 
compared with experimental measurements where they are 
known. All values are expressed in %. 



Parent 2002 
nucleus value [Jl] 



This work 



unsealed scaled Adopted 



expt 





= 0: 










38 J^m 


0.090(30) 


0.068 


0.062 


0.085(30) 


< 0.28" 


^^Sc 


0.020(20) 


0.007 


0.027 


0.015(15) 


0.040(9)' 


46y 


0.035(15) 


0.008 


0.024 


0.025(20) 


0.053(5)" 




0.045(20) 


0.049 


0.047 


0.040(20) 


< 0.016" 


s-^Co 


0.040(20) 


0.049 


0.038 


0.050(20) 


0.035(5)" 


e^Ga 


0.085(20) 


0.160 


0.159 


0.120(40) 


< 0.040(15)" 


66 As 


0.020(20) 


0.110 


0.087 


0.050(30) 




™Br 


0.070(20) 


0.226 


0.235 


0.150(80) 
0.050(30)'' 






0.050(30) 


0.045 


0.047 


< 0.075" 



"From Hagberg et al. (1994) [SJ 

''From Daehnick and Rosa (1985) [ssll averaged with earlier re- 
sults. 

<^From Hyland et al. (2006) O 

''No new calculations were performed for ^''Rb. 

"^From Picchaczck et al. (2003) [H 



isfactory. But in the upper p/-shell, the calculated value 
for ^^Ga is three times larger than measured in recent 
experiments [l5|. Shell- model calculations in this region 
are complicated by the massive size of the Hamiltonian 
matrices. To keep our calculations tractible, we kept the 
/7/2 shell closed in these cases, but there is considerable 
evidence [29| that this could be a poor assumption. 



IV. THE RADIATIVE CORRECTION 



Prior to 1990 



Conventionally, the radiative correction has been sep- 
arated into two parts, one that contains the nucleus- 
dependent terms, called the 'outer' radiative correction, 
and one that is independent of the nucleus, the 'inner' 
radiative correction. Principally due to the work of Mar- 
ciano and Sirlin (for example, Refs. [37l. [38l. [39| ) . the ra- 
diative correction applied to the uncorrected /3-decay rate 
was expressed as follows: 



p = r°(i + <5y(i + A^) 



(23) 



^ In a few cases, the state calculated to have the largest charge- 
dependent admixture was the second excited 0+ state. In these 
cases we optimized the agreement between theory and experi- 
ment for the excitation energy of that state 



a 
2^ 
a 
2^ 



[9iErn)+S2+S3 

( mp 



3 In 



V 2i?„ 



large E„ 



81 47r2 ^ 

h (^2 + <53 

10 3 



(24) 
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a;: 



2^ 



OL 

2^ 



31n- 



41n- 



■+ln 
-4 In 
+ ln- 



m, 
m 

TO 

TO. 



w 



TO„ 



f 2C 
+ 2C + A 



(25) 
(26) 



where Em is the maximum electron energy in /3-dccay, 
and TO„, TOp, TOz arc the masses of the W^-boson. proton 
and Z-boson. The separation into outer and inner terms 
is accommodated in 8'^ and AJJ^ respectively. 

In the outer correction. (5^, the order-a term contains 
the function 'g{Em)- it is the average over the beta en- 
ergy spectrum of the function g{E, Em.Y which was de- 
fined by Sirlin (see Eq. (20b) of Ref. [331) and is not re- 
produced here. Its large- limit is shown in Eq. ([M]), 
indicating that the expression is dominated by the log- 
arithm, ln(TOp/(2i?„i)). The last two terms in the outer 
correction. 82 and ^3, represent corrections to order Za^ 
and Z^a'^ respectively. The origin of the g[E,n) term 
- together with that of the leading term in the inner 
radiative correction, 31n(TOw/TOp) - is the jW-hoic and 
bremsstrahlung diagrams, which are taken together to 
remove the divergence as the photon energy goes to zero. 
Both $2 and ^3 also come from a standard QED calcula- 
tion of the 7Vl^-box and bremsstrahlung graphs [iO, [4l[ , 
but in their case the electron was allowed to interact with 
the Coulomb field of the nucleus. Care was taken not to 
double count with the Fermi function. The calculation 
was complete to order Za^ but only estimated in order 

In the inner correction, A^, the second and third 
terms, ln(TOw/TOA)+2C, like the first term, also represent 
a 714^ box graph, but this time it involves an axial- vector 
weak interaction. The evaluation of this graph can be di- 
vided into two energy regimes: the high-energy (or short- 
distance) part given by the logarithm, and the low-energy 
(or long-distance) part denoted by 2C. The parameter 
TOa, referred to as the low-energy cut-off, divides these 
two energy regimes. Marciano and Sirlin [sll allowed it 
to take on a range of valueSj400 MeV < TOa < 1600 MeV 
(revised slightly by Sirlin [44I to be TOai/2 < TOa < 2TOaj, 



with 



being the Ai-vector-meson mass). The low- 



energy component, 
contribution 



2C, was approximated by its Born 



C ^ CBorn = 3.gA (0.266) (mp + f^n) = 0.885, (27) 

where gj^ = 1.26 is the axial vector coupling constant 
accepted at the time and {fip + fin) = 0.88 is the nuclcon 
isoscalar magnetic moment. The factor 0.266 is the value 
of the loop integral that was rendered finite by the use 
of dipole form factors for the nucleon electromagnetic, 
7A^, and axial-vector, WN, vertices. The fourth term 
in Eq. (|25[) . with the logarithm ln(TOw/TOz), arises from 
ZW-hoyi graphs; while the last term, Ag, represents a 
small perturbative QCD correction that was evaluated 
by Marciano and Sirlin [s^ to be An = —0.34. 



The value of the outer radiative correction as defined 
in Eq. ((M]) . ranges from 1.39-1.65% for the known super- 
allowed emitters (see Ref. Ql)- Following Sirhn [i^l, the 
assigned uncertainties are set equal to {a/2T:)5z as an es- 
timate of the error made in stopping the calculation at 
that order. The value of the inner radiative correction as 
obtained from Eq. ^ with C from Eq. ^ is [H, [13] 



(28) 



K{old) = 2.40(8)%. 



These results provide the essential foundation of the 
radiation corrections still used today. However a number 
of improvements have been introduced in the intervening 
17 years. 



B. A nuclear-structure dependent term 

The low-energy part of the 7Vl^-box diagram for an 
axial-vector weak interaction, denoted 2C, was approxi- 
mated by its Born contribution in Eq. (j27p , and was eval- 
uated on a single nucleon. However, in a finite nucleus 
with many nucleons present, Jaus and Rasche (43j ob- 
served that the two hadronic-interaction vertices, 7iV and 
WN , do not have to be with the same nucleon. Thus, 
in finite nuclei there can be two types of contributions: 
those in which 77V and WN vertices are with the same 
nucleon and those in which they are not. The evaluation 
of the former terms yields expressions |44| that are pro- 
portional to r+, the isospin ladder operator, and so are 
also proportional to the Fermi /?-decay operator. There- 
fore, they produce a universal correction - the same in all 
nuclei - with the value Cboih, which is given in Eq. (|27p . 
The remaining terms, those in which the interactions 
are with different nucleons, must be evaluated with two- 
body operators that depend on the nuclear structure of 
the states involved. Thus, the expression for C given in 
Eq. (|?7|) must be replaced by the following equation: 



C = C, 



Born 



c. 



NS, 



(29) 



where Cns comprises the nuclear-structure dependent 
terms. Calculations of Cns were first made in 1992 [3, 

m- ^ 

A further modification was introduced in 1994 [46[. 
In calculations of Cboih that had been made up to that 
time, the axial-vector and electromagnetic coupling con- 
stants, (7a and (fj.p + ^„) - see Eq. (j27p - had been given 
their free-nucleon values. Yet there is ample evidence 
in nuclear physics that coupling constants for spin-flip 
processes are quenched in the nuclear medium, with 
the amount of quenching varying from nucleus to nu- 
cleus. Thus, one should really be replacing C^^^^, the 
value obtained with free-nucleon coupling constants, with 



ucnchcd 
Born 



However, to separate the nucleus-dependent 



and nucleus-independent parts of the latter, we write 



quenched 
Born 



^*^Born 



free 
Born 



+ iq- m 



free 



(30) 
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where q is the factor by which the product of the weak 
and electromagnetic coupUng constants is reduced in the 
medium relative to its frec-nucleon value. 

The first term in Eq. (|30p . which remains universal, is 
retained in the inner radiative correction, replacing C 
in Eq. (|25p . The second term becomes part of a sep- 
arate nuclear-structure-dependent radiative correction, 
SnS: which also includes C^^g"''^'''^, the value of Cns re- 
calculated with quenched operators. This correction is 
written as 



'JNS 



^qu^chcd ^ _ ^^^fi 



(31) 



and is incorporated with the other nuclear-structure- 
dcpendent correction term, Sc - see Eq. ([S|). Calcu- 
lated values of Sns range from -0.360% to +0.030%, 
each generally being smaller in magnitude than the cor- 
responding value of Sc- 

We return to 6ns in Sect. lIVEl 



where NLL is a next-to-leading log correction that Czar- 
necki et al. estimate to be NLL = —0.0001. The coeffi- 
cient a(m) is a running QED coupling constant whose 
value at m = jrip is 1/133.986 and at to = TOe is 
1/137.089 [13. 

This new result can still be organized to pre- 
serve the separation of nucleus-dependent and nucleus- 
independent components. The separation we hereby 
adopt is 



1 + 5'r 



a 



2Er, 



[L{2Era,mp) + ^[52 + h]] (37) 



5(to„to,) + -<- 



2tt 



TT 

TO, 



NLL. (38) 



C. Improvements to 

In 2005, Czarnecki. Marciano and Sirlin [4^ revisited 
the 0(a^) correction for neutron beta decay. They began 
by trivially updating the value of C^^°^ to reflect the 
current value of the axial- vector coupling constant, = 
1.27, to get 



free 
Born 



0.891, 



(32) 



which replaces the value given in Eq. (P7)) . 

They then went on to re-evaluate A^, focusing par- 
ticularly on the leading log corrections. U sing an es- 
tablished renormalization group summation [39| for the 
leading short-distance logs, S{'mp,mz), they extended 
the method to the lower energy region between 2Em and 
nip to obtain L{2E„i,mp). This resulted in the replace- 
ments 



1 



— In — 

TT TO.„ 



ia , TO, 
— In — - 

2-K 2E, 



S'(TOp,TOz) = 1.02248 (33) 
L{2E^,mp), (34) 



where 



L{2Em,mp) = 1.026725 



2a(TOe), 2E„ 
m 



-1 9/4 



37r 



(35) 

The complete radiative correction, RC, including order 
Za^ and Z'^a^ terms, could then be written 14711 



1 + RC 



1 + ^ 
2tt 



31n- 



2E„ 



{l(2£;„. Top) + £ [2Cf5-^„ + 62 + 6,] } 
a{mp) 



X <^ S{mp,mz) + 



27r 



1 "^D A 

In -^+Aq 


+ NLL^ 







(36) 



We will use this separation here and in our future work 
on superallowed (5 decay. It results in a small change to 
the values of (5^ and that we used in our recent review 



D. Reduced uncertainty for Ap 



In Sect. lIV.^ we explained that the terms 
ln(TOw/mA) + 2C in Eq. p5)) arose from the ^W- 
box graph for an axial-vector weak interaction. These 
two terms came from splitting the evaluation of this 
graph into two energy regimes. The division between 
the two regimes was chosen to be to a = 1.2 GeV (4^ . 
roughly the mass of the Ai resonance, and its range 
of uncertainty was taken to be from TOa/2 to 2toa. 
This ad hoc range determination actually produced the 
largest single contributor to the uncertainty in the CKM 
matrix element, Vud- 

To reduce the hadronic uncertainty in the radiative 
correction, Marciano and Sirlin Q have looked again at 
the 7VF-box graph for an axial-vector weak interaction. 
This time they split it into three energy regimes, rather 
than two and, where possible, they drew on independent 
information to control their results: 

• Short distances, (1.5 GeV)^ < < oo: This is 
a domain where QCD corrections remain pertur- 
bative. Marciano and Sirlin added higher-order 
terms, noting that these terms are identical (in the 
chiral limit) to QCD corrections to the Bjorken sum 
rule for polarized electroproduction and can there- 
fore be obtained from well-studied calculations for 
that process. 

• Intermediate distances, (0.823 GeV)^ < < 
(1.5 GeV)^: In this region, they used an inter- 
polation function between low and high energies, 
motivated by vector-meson and axial- vector-meson 
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TABLE V: Calculated transition-dependent radiative correc- 
tion, 5'fi, in percent units, and the component contributions. 
In our previous works {eg. Ref. [j) 5^ was defined as the sum 
of the contents of columns 2-4; this result is given in column 
5 and labeled "Former (5^." As explained in the text, we have 
now redefined 6'^ to include the additional term in column 6; 
the new values for 5'^ are given in the last column. 



Parent 

mucleus ^^(Sr? 



Former 
5'r 



Redefined 



T — 


1 . 

— ± . 






















10^' 


1 

± , 








n 

u 




1 

_L 


fi^9 


n 

u. 


097 




14 Q 


1 

± , 


. ^ou 


n 

u 


99fi 


n 

u 




1 

J_ 


"190 


n 

u. 




1 ^4'^CS"! 


IN c 


1 

± , 




n 


9fi8 


n 

u 


ni 9 


1 

_L 


.4:04: 


n 


099 




iVlg 


1 

± , 


1 99 


n 

u 


.OU 1 


n 

u 


ni 7 


1 

_L 


44fi 


n 

u. 


n9n 




26c: 
Ol 


1 

1 , 






■^49 

.04:^ 


n 

u 


09"^ 
.uzo 


1 

J_ 


49n 


n 

u. 


ni Q 




30 g 


1 


.005 


0, 


.371 


0, 


.029 


1 


.405 


0, 


.018 


1.423(29) 


3^Ar 


0, 


.963 


0, 


.396 


0, 


.035 


1 


.395 


0, 


.017 


1.412(35) 


^«Ca 


0, 


.929 


0, 


.426 


0, 


.042 


1 


.397 


0. 


.017 


1.414(42) 


42^. 


0, 


.906 


0, 


.456 


0, 


.050 


1 


.412 


0. 


.016 


1.428(50) 


T, = 


0: 






















26 


1, 


.110 


0, 


.328 


0, 


.020 


1 


.458 


0, 


.020 


1.478(20) 


3*01 


1, 


.002 


0, 


.390 


0, 


.032 


1 


.425 


0. 


.018 


1.443(32) 


38 771 


0, 


.964 


0, 


.420 


0, 


.039 


1 


.423 


0, 


.017 


1.440(39) 


^^Sc 


0, 


.939 


0, 


.451 


0, 


.047 


1 


.436 


0, 


.017 


1.453(47) 


46y 


0, 


.903 


0, 


.472 


0, 


.054 


1 


.429 


0, 


.016 


1.445(54) 


50Mn 


0, 


.873 


0, 


.494 


0, 


.062 


1 


.430 


0, 


.015 


1.445(62) 


s^Co 


0, 


.844 


0, 


.513 


0, 


.071 


1 


.428 


0. 


.015 


1.443(71) 


62 Ga 


0, 


.805 


0, 


.553 


0, 


.087 


1 


.445 


0. 


.014 


1.459(87) 


66 As 


0, 


.791 


0, 


.570 


0, 


.095 


1 


.456 


0. 


.014 


1.470(95) 


™Br 


0, 


.776 


0, 


.591 


0, 


.105 


1 


.473 


0, 


.013 


1.49(11) 


^%b 


0, 


.761 


0, 


.609 


0, 


.115 


1 


.485 


0. 


.013 


1.50(12) 



dominance. By limiting the number of terms to 
three, they had sufficient matching conditions to 
determine the coefficients uniquely. 

Long distances, < < (0.823 GeV)^: 

Integrating the long-distance ampUtudes up to 
= (0.823 GeV)^, where the integrand matches 
smoothly to the interpolation function, they ob- 
tained a smaller value for CS°° , 



free 
Born 



0.829, 



(39) 



than given in Eq. ([32]) . However this smaller value 
is caused entirely by the reduction in the effective 
upper limit to the loop integration, and is almost 
completely compensated for by the consequently 
higher values obtained for the graph in the other 
energy regimes. 

In the end, Marciano and Sirlin Q find that the net 
effect of this re-evaluation of the 7W^-box axial graph 
is a very small reduction in the radiative correction of 
1.4 X 10""'. More important than this reduction, the 
new method provides a more systematic estimate of the 



hadronic uncertainties. Allowing for a ±10% uncertainty 
for the CgQ^j correction in Eq. (|39p . a ±100% uncertainty 
for the interpolator contribution in the intermediate re- 
gion, and ±0.0001 uncertainty from neglected higher or- 
der effects, Marciano and Sirlin Q find the total uncer- 
tainty in the radiative correction is ±0.00038. This cor- 
responds to more than a factor of two reduction in the 
loop uncertainty for hadronic effects (c/. Eq. ipS]) ). 



E. New values for S'^, and 5ms 

We maintain the traditional separation of the radiative 
correction into a nucleus-dependent outer correction and 
a nucleus-independent inner correction - see Eqs. (j37p 
and ([55]) . This means that the outer correction, (5^, is 
slightly redefined and is now written as 



27r 



[g{E,n) + 62 + 63 + 6^ 



(40) 



where the new term, , simply represents the differ- 
ence between the definition of S'j^ given in Eq. (|37p and 
that given in Eq. p4|) . It is the leading- log extrapolation 

which is contained in the 



of the logarithm hi{mp/2Em) 
function g{Ejn)- Values of Sa 



and the redefined are 



TABLE VI: Calculated nuclear-structure-dependent radiative 
correction, Sns, in percent units, and the component contri- 
butions. 





2002 


This work 


Parent 


Sns{%) 


Sns{%) 


nucleus 


Ref. \£ 


^.^ucnched (g_i)c7|,- adopted 





-1: 










-0.360(35) 


-1.318 


-0.176 


-0.345(35) 


14q 


-0.250(50) 


-0.844 


-0.208 


-0.245(50) 


i«Ne 


-0.290(35) 


-1.051 


-0.198 


-0.290(35) 


22 Mg 


-0.240(20) 


-0.750 


-0.213 


-0.225(20) 


26Si 


-0.230(20) 


-0.705 


-0.227 


-0.215(20) 


30g 


-0.190(15) 


-0.557 


-0.242 


-0.185(15) 


34 Ar 


-0.185(15) 


-0.520 


-0.257 


-0.180(15) 


3«Ca 


-0.180(15) 


-0.475 


-0.271 


-0.175(15) 


42^i 


-0.240(20) 


-0.765 


-0.241 


-0.235(20) 


= 


0: 

0.009(20) 


0.242 


-0.227 


0.005(20) 


^^Cl 


-0.085(15) 


-0.118 


-0.257 


-0.085(15) 


38 j^m 


-0.100(15) 


-0.158 


-0.271 


-0.100(15) 


42Sc 


0.030(20) 


0.391 


-0.241 


0.035(20) 


46y 


-0.040(7) 


0.093 


-0.248 


-0.035(10) 


50Mn 


-0.042(7) 


0.084 


-0.254 


-0.040(10) 


«*Co 


-0.029(7) 


0.112 


-0.261 


-0.035(10) 


62Ga 


-0.040(20) 


0.087 


-0.272 


-0.045(20) 


66 As 


-0.050(20) 


0.010 


-0.278 


-0.060(20) 


^OBr 


-0.060(20) 


-0.085 


-0.283 


-0.085(25) 




-0.065(20) 


-0.026 


-0.288 


-0.075(30) 
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TABLE VII: Corrected J^t values for the thirteen best known 
superallowed decays, obtained with the new correction terms 
presented in this work. The experimental ft values were taken 
from results in our 2005 survey [j| updated with more recent 
published data 0, Q, i, [13, 0, IllH [H, El, Eil • The average 
Tt value and the normalized of the fit to a constant appears 
at the bottom. 



Parent 






nucleus /i(s) 


5'n{%) Sns{%) 5c{%) 





= -1: 

^"C 3039.5(47) 1.679(4) 

3042.5(27) 1.543(8) 

^^Mg 3052.2(72) 1.466(17) 

Ar 3052.5(82) 1.412(35) 



34 



-0.345(35) 0.175(18) 3074.5(49) 

-0.245(50) 0.330(25) 3071.6(33) 

-0.225(20) 0.380(22) 3078.3(74) 

-0.180(15) 0.665(56) 3069.4(85) 



3^01 

38 J^T72 

^^Sc 

46y 

s^Mn 
5*Co 
^^Ga 
Rb 



74 



= 0: 

3037.0(11) 
3050.0(11) 
3051.1(10) 
3046.4(14) 
3049.6(16) 
3044.4(12) 
3047.6(15) 
3075.5(14) 
3084.3(80) 



1.478(20) 
1.443(32) 
1.440(39) 
1.453(47) 
1.445(54) 
1.445(62) 
1.443(71) 
1.459(87) 
1.50(12) 



0.005(20) 
-0.085(15) 
-0.100(15) 
0.035(20) 
-0.035(10) 
-0.040(10) 
-0.035(10) 
-0.045(20) 
-0.075(30) 



0.310(18) 
0.650(46) 
0.655(59) 
0.665(56) 
0.620(63) 
0.655(54) 
0.770(67) 
1.48(21) 
1.63(31) 



3072.5(15) 
3071.3(21) 
3071.7(24) 
3071.2(27) 
3073.4(30) 
3066.9(28) 
3066.7(33) 
3073.0(72) 
3077(13) 



Average J^t 

2 / 



3071.4(8) 
0.6 



given in Tabic |V] for all the superallowed transitions of 
interest. 

The new inner correction is defined by Eq. (|38|) . with 
^Born taken from Eq. p9|) . With its uncertainty obtained 
from Marciano and Sirlin H , the result is 



A;^ = (2.361 ± 0.038)%. 



(41) 



It is important to note that with the re-evaluation 
of Cg°°jj, there is a consequent change in the nuclear- 
structure dependent correction 5ns given in Eq. 
Fortunately, the change is very small, being 



new /--^old 
Born ^ "-"l 



Born)(":) 



-0.3(-0.062)2.3 X lO'^' 
0.004%. (42) 



In addition to making this change, we have also taken the 
opportunity to re-evaluate Cns using the more recently 
available shell-model effective interactions described in 
Sect, nil A 31 Our revised 5ms values are listed in Ta- 
ble IvH As in Tablesim and lllli we give (in columns 3 and 
4) the results from one sample calculation for each nu- 
cleus. Then in column 5 we present adopted 5ns values 
that result from our assesment of all calculations made 
for each decay, not just the ones shown in columns 3 and 
4; the uncertainties were chosen to encompass the spread 
in the results from those calculations. For comparison, in 



column 2 we list the values we adopted for S^s in 2002 
[1]. In all cases the new values agree with the old ones 
within the quoted uncertainties. 



V. Tt VALUES, Vud AND CKM UNITARITY 

We have calculated improved results for the correction 
terms 5ci (see TableEII, 5c2 (Table|Tl| and 5ms (Ta- 
ble l Vip : and, based on the work of Marciano and Sirlin, 
we have presented revised values for 5'j^ (Table|V| and 
(Eq. (|4ip ). We are now in a position to extract corrected 
J-t values from the current world data for superallowed 
0+ 0"*" transitions. 

We use the same data set as that described in Sect.llll 
it represents an interim update of our 2005 complete sur- 
vey PI and includes ten additional published measure- 
ments 0, 0, i, [13, [m, E m, H E [ll- Results are 
given in Table lVTIl for the thirteen superallowed transi- 
tions whose ft values are known to a precision of 0.3% 
or better. The Tt values given in column 6 were ob- 
tained from the data in the preceeding columns through 
the application of Eq. The corrected Tt values are 
also plotted in Figure [2] 

It is clear from the normalized given on the bottom 
line of the table that the statistical agreement among the 
J-t values remains excellent. Furthermore, it is evident 
from the figure that no longer shows any deviation 
from the overall average as it did in Fig.[TJ However, it is 
equally evident that instead the ^°Mn and ^"'Co Tt values 
are now low, and by amounts that are no less statistically 
significant than the amount by which the "^^V value was 
previously high. 

Rather than being a negative result, however, this pos- 
sible discrepancy offers us the opportunity to use the 
cases of ^"Mn and ^''Co as a valuable test of our im- 
proved calculations. The Qec value for each of them 



3100 
3090 
'in 3080 

3070 
3060 



"O "Al" "CI «'Ga 
"Ar "Sc "Mn"Co 



5 10 15 20 25 30 35 

Z of daughter 

FIG. 2: Results for the new corrected Tt values (from Ta- 
ble lVII|) for the thirteen best known superallowed decays. The 
corresponding uncorrected ft values appear in the left panel 
of Fig.[T] The shaded horizontal band gives one standard de- 
viation around the average Tt value. 
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has been measured only twice with (claimed) high pre- 
cision [4^, [4^ , and one of these references [31 also in- 
cluded a measurement of the Qec value for ''^V, which 
Pcnning-trap measurements have recently shown [ll| 
to be low by 2 keV - more than three times its originally 
quoted standard deviation. If, as seems likely, the prob- 
lem with the measurement in Ref. [i^ is not limited 
to that measurement alone, then doubt is certainly cast 
on the ^°Mn and ^^Co Q^c-value results quoted in that 
reference as well. New Penning-trap measurements of 
both Qec values are currently in progress (HO], and the 
question should be settled shortly. If the Qec values in 
Ref. [4^ prove to have been too low again, then the new 
Penning-trap measurements will serve to increase the Tt 
values for ^°Mn and ^^Co and could well bring them into 
close agreement with the average Tt value. If so, this 
would add strong support to our new calculations. 

The average corrected Tt value obtained from our new 
analysis, 3071.4(8) s, is lower by more than one standard 
deviation, compared to the comparable result obtained 
in our 2005 survey, 3072.7(8) s. If the new measure- 
ments do prove to increase the Qec values for ^°Mn and 
^^Co, then this discrepancy will decrease slightly, but 
there is no avoiding the fact that the inclusion of some 
core orbitals in the nuclear-structure-dependent correc- 
tion terms has increased the correction in a number of 
cases, which in turn leads to a reduction in their Tt val- 
ues. A significant change in the nuclear model has led to 
a significant change - but not a revolutionary one - in 
the average Tt value. 

The new average J-t value yields a new value for Vud 
via the equation 



2Gl{l + Al)Tt' ' ' 

where Gp is the well known weak-interaction constant 
for the purely leptonic muon decay [l7j . It has been our 
practice when using the Tt value in this context to add 
0.85(85) s to its value to account for possible systematic 
errors in the treatment of the radial wave function in the 
calculation of 8c- (This point is discussed in detail in 
section III C of Ref. [if.) Continuing this practice, we 
obtain the following result for the up-down element of 
the CKM matrix: 

IKdl = 0.97418(26). (44) 

This result can be compared with the value 0.97380(40), 
which was obtained in 2005 [l[. The new value is (just) 
within the uncertainty of the previous value, and carries 
an uncertainty that is one third smaller. 

The final step is to combine this new value of iKidl with 
the other top-row elements of the CKM matrix, \yus \ and 
\yuh\i to test the unitarity of the matrix. Taking the 
values of the latter two elements from the 2006 Particle 
Data Group review we obtain the stunning result 

|Kd|' + |K.|' + |Kfa|' = 1.0000 ±0.0011. (45) 



Unitarity is fully satisfied with a precision of 0.1%. 



VI. CONCLUSIONS 

We have presented new calculations of the nuclear- 
structure-dependent corrections to superallowed 0^ — *■ 
0+ nuclear /3 decay. The calculations incorporate core 
orbitals in the shell model in cases where independent ex- 
perimental information indicates that they are required. 
Where possible, they also make use of effective interac- 
tions that have been published since our previous cal- 
culation of these correction terms As in that work, 
we have included twenty transitions in our calculations, 
thirteen that are by now rather well measured and seven 
more that are likely to be accessible to precise measure- 
ments in the future. 

The agreement among the corrected Tt values for the 
thirteen well measured cases is very good, although there 
is a possible small discrepancy for the cases of ^°Mn and 
^^Co. A new Penning-trap measurement of the Qec val- 
ues for these two transitions is expected in the near fu- 
ture, and its effect on this discrepancy could serve to test 
the validity of our calculations. 

With our new corrections, the value of |V^d| is in- 
creased by 0.04%, or by one standard deviation of the 
previous result hi. With the new value, the sum of 
squares of the top-row elements of the CKM matrix is 
in perfect agreement with unitarity. 

The improved calculations presented in this work were 
inspired by the remarkable recent improvements in exper- 
imental precision, particularly in the measurement of the 
"'^V Qec value. The only way that the calculated cor- 
rections can be tested and improved is by such precise 
measurements, both on the currently well-known transi- 
tions and on other as-yet-unstudied superallowed transi- 
tions that have larger calculated corrections. If the cal- 
culated correction terms replace the significant scatter in 
the measured /i values (see the left panel in Fig.[T|) with 
a set of self-consistent corrected Tt values, then they can 
surely be relied upon to produce a secure value for |U„rf|. 
The present calculations testify to the value of increased 
experimental precision. 
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